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Abstract 

Previous work has shown that the standard supergravity approximation can break down when 
using AdS/CFT duality to study certain top-down formulations of the jet stopping problem in 
strongly-coupled AA=4 super- Yang-Mills (SYM) plasmas, depending on the virtuality of the source 
of the "jet." In this paper, we identify the nature of this breakdown: High-momentum gravitons 
in the gravitational dual get stretched into relatively large classical string loops by tidal forces 
associated with the black brane. These stringy excitations of the graviton are not contained in 
the supergravity approximation, but we show that the jet stopping problem can nonetheless still 
be solved by drawing on various string-theory methods (the eikonal approximation, the Penrose 
limit, string quantization in pp-wave backgrounds) to obtain a probability distribution for the 
late-time classical string loops. In extreme cases, we find that the gravitons are stretched into very 
long folded strings which are qualitatively similar to the folded classical strings originally used by 
Gubser, Gulotta, Pufu and Rocha to model the jet stopping problem. This makes a connection 
in certain cases between the different methods that have been used to study jet stopping with 
AdS / CFT and gives a specific example of a precise AA=4 SYM problem that generates such strings 
in the gravity description. 



I. INTRODUCTION 



A. Background 

Inspired by the observation of (and rapidly growing body of experimental information on) 
jet quenching in relativistic heavy ion collisions, there has for many years been an interest 
in the theory of jet quenching and what can be learned about that theory by studying 
interesting limiting cases. One of the simplest-to-pose thought experiments is this: How 
far does a very-high momentum excitation (the potential precursor of a would-be jet) travel 
in a thermal QCD medium before it loses energy, stops, and thermalizes in the medium? 
And how does the answer to that question depend on the effective strength a s of the strong 
coupling? 

This question can be addressed from first principles in various theoretical limits. One such 
limit is that of weak coupling, which in principle applies to asymptotically large temperatures 
T and jet energies E, for which the relevant running values of a s are small. In that limit, 
the stopping distance £ stop for a high-energy parton (E S> T) scales with energy as E 1 / 2 , up 
to logarithms. 1 A contrasting limit of interest occurs when the running values of a s relevant 
to jet stopping are all large. 2 This problem is not very tractable from first principles in 
QCD itself, but, through gauge-gravity duality, progress can be made for QCD-like plasmas 
with gravity duals, such as A/=4 super Yang Mills (SYM) theory. For some years, people 
have considered various ways to study analogs of jet stopping in such plasmas, namely the 
stopping distance for various types of localized, high-momentum excitations. The exact 
stopping distance depends on details of exactly how the "jet" is prepared, but universally 
these studies have found that the maximum possible stopping distance £ max scales with 
energy as E 1 / 3 0-10], in contrast to the weak-coupling scaling of E x l 2 . This is an interesting 



theoretical result because it teaches us that the scaling of jet stopping with energy depends 
on the strength of the coupling. It remains an open question (which we will not answer 
here) how E 1 / 3 starts to move toward E 1 ^ 2 as one lowers the coupling, and vice versa. 3 

The stopping distance of high-momentum, localized excitations traveling through the 
plasma depends on more than just the energy of the excitation. Depending on exactly how 
one creates the excitation (the "jet"), one may get stopping distances £ stop significantly 
smaller than the maximum £ max . As an example from weak coupling, imagine that we 
spread out the total energy and momentum E of the jet among 10 partons, each having 
energy E/10, rather than putting it all into a single parton of energy E. Each of the 10 
partons has lower energy than the single one and so will stop sooner; so the stopping distance 



1 A specific weak-coupling calculation of the stopping distance for QCD in the high-energy limit may be 
found in ref. p|. However, the scaling of this result was implicit in the early pioneering work of refs. [2, S] 
on bremsstrahlung and energy loss rates in QCD plasmas. Introducing supersymmetry will not change 
the conclusion that the stopping distance scales as E 1 / 2 (up to logarithms) at weak coupling. 

2 In a weak-coupling analysis, the two running couplings relevant to jet stopping are, roughly, a s (T) and 
a s (Qj_), where Q±_ ~ (qE) 1 / 4 grows slowly with energy and is the scale of the typical relative momentum 
between two daughter partons when a high energy parton splits through hard bremsstrahlung or pair 
production, (q ~ ot 2 T^ is a scale characteristic of the plasma that parametrizes transverse momentum 
diffusion of high-energy partons.) A third limiting case of interest, not addressed in this paper, is where 

a s (T) is large but a s (Q±) isLsmall. See, for example, Liu, Rajagopal, and Wiedemann 
See the conclusion of ref. for further discussion of this point. 
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FIG. 1: Examples of the decay of a very high-energy (a) slightly virtual photon, (b) slightly virtual 
graviton, or (c) on-shell W + boson, inside a standard-model quark-gluon plasma, producing high- 
momentum partons moving to the right. In the context of W=4 super Yang Mills, the q, u, and d 
above represent adjoint-color fermions or scalars carrying R charge. For strong coupling, of course, 
one should not picture perturbatively, as in this figure, the high-momentum excitation created in 
the plasma by the decay. 



for the high-momentum excitation depends on how many high-energy partons we use in the 
initial state. In the weak-coupling case, the maximum stopping distance £ max corresponds to 
the particular initial state where all the energy is concentrated into a single initial parton. 

In the strong-coupling case, we cannot speak of individual partons, but the stopping 
distance a gain depends on how we prepare the initial high-momentum excitation. In our 
work (l. [9I. . we create the initial excitation in a way that is analogous to what you would 
get if a high-momentum, slightly virtual photon (or graviton or other massless particle) 
decayed hadronically in the quark-gluon plasma, as depicted in fig. [TJ Alternatively, one 
could consider the decay of a high- momentum on-shell W boson (also depicted). For these 
methods of creating "jets," one finds that the maximum possible stopping distance scales as 

£1/3 

C ax ~^p- (1.1) 

As we will review later, it turns out that the stopping distance may be made smaller than 
(11. ip by varying the virtuality — q 2 = —Qn(f_ of the virtual photon (or equivalently the mass- 
squared M 2 of the on-shell W boson) [g, |9|. The important point is that there is a range 
of stopping distances 4to P ^ Cax for our "jets," depending on the details of how those 
excitations are created. 

Most top-down studies of jet stopping using gauge-gravity duality have studied the infi- 
nite color and infinite coupling limit, N c =oo and A=oo, where A = N c g\ M is the 't Hooft 
coupling. To understand the true high-energy behavior, however, it is important to study 
the corrections to these limits. As an example, fig. [2] shows two different scenarios one might 
imagine for the maximum stopping distance £ max for strongly-coupled A/"=4 SYM. One is 
that £ max grows like £ 1/3 at high energy, up to arbitrarily high energies. The other is that 
it starts growing like £ 1/3 at high energy E 3> T, but then crosses over to some different 
power-law behavior once E exceeds some positive power of A times T (e.g. A 2 T in the figure). 
In the latter case, E 1 ^ 3 would not be the true behavior for arbitrarily large E and large but 
finite A. But there is no way to tell the difference between these two scenarios if one only 
has A=oo calculations! For this reason, three of us analyzed the parametric size of finite-A 
corrections to jet stopping distances in ref. [H| . We found that the formal expansion in 
1 / vA (which corresponds to an expansion in the string parameter a' on the gravity side) 
breaks down for some jets and is safe for others, depending on the stopping distance 4top 
of the jet (and therefore on the virtuality — q 2 ), as depicted in fig. [3j The stopping distance 
and virtuality parametrize the horizontal axis in this figure. The vertical axis is a relative 
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FIG. 2: Examples of two different scenarios for the high-energy (E ^> T) behavior of the maximum 
jet stopping distance £ max (E) which are indistinguishable with A=oo calculations. 



measure of the importance of a given correction compared to the A=oo result (see ref. [11 
for details). The curves are labeled by the sequence of higher- curvature terms in the gravita- 
tional dual theory action that correspond, via the AdS/CFT correspondence, to a sequence 
of corrections in powers of 1/vA in the 3+1 dimensional J\f=4 SYM theory. Throughout, 
N c is taken to be infinite. The result of this study was that, for A > 1, corrections to the 
A=oo result are parametrically small for A~ 1/,6 £ max <C £ st0 p IS ^max- In particular, corrections 
to the maximum stopping distance £ max ex E 1 ^ are small. But the interesting case is when 
jets are created in such a way that 

T- 1 < 4to P < A~ 1/6 £ max , (1.2a) 

which is ^ g 

T- 1 « 4to P < ^(ff ■ (1.2b) 

In this case, the fate of A=oo results for the stopping distance was unclear. For £ stop ~ 
A _1 / 6 £ max , all the corrections are the same size, and so the formal expansion in powers of 
1 / \f\ has broken down. Yet the individual corrections are all small (of relative importance 
A _1//2 ) for that 4top- From fig. [3j we cannot tell whether the sum of the corrections to A=oo 
will remain small for £ stop < A _1 / 6 £ max or whether, instead, the A=oo calculation becomes 
useless there. 

The purpose of the present paper is to understand the physics (on the gravity side) of 
what is going on in the region ( II. 2p where the naive expansion in powers of 1/ \f\ (powers 
of a') breaks down, and to figure out how to account for the effect of this physics on the 
jet stopping distance. Note that the interesting window (1 1 . 2a[) of stopping distances exists 
only if the energy is large enough that T _1 <C A _1 / 6 £ max . By ( II. ip . this requires E ^> A X / 2 T, 
which we will assume throughout the rest of this paper. 

Before outlining what we have done, it will be useful to first explain one other quali- 
tative feature of the A=oo calculation. Excitations created in the field theory correspond 
to excitations created on the boundary of AdSs-Schwarzschild, which then fall towards the 
black brane over time, such as depicted in fig. HI The 3-space distance that this excitation 
travels before falling into the horizon matches the stopping distance of the corresponding 
excitation in Af=4 SYM. 4 For £ stop 4nax, which includes the region (II. 2p of interest, there 



4 See ref. [9( for a discussion in the context of the present paper, but this correspondence is implicit in the 
earlier work of refs. 

00. 
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-q 2 /T 2 < 1 1 1— 

X p/3 A 2/3£.2/3 £.2/3 



FIG. 3: A parametric picture of the relative importance of higher-derivative corrections to the 
low-energy supergravity action as a function of the stopping distance ^stop 

(using the A=oo result 

for £ s top)- The axis are both logarithmic, and an importance of 1 indicates that the individual 
correction would, by itself, significantly modify the A=oo analysis. The measure of "importance" 
is explained in ref. Also shown, as an alternative horizontal axis, is the 4-dimensional virtuality 
— q 2 of the source that created the jet, where E = E/T. 



is a nice simplification. On the gravity side, the excitation falling in fig. H] turns out to 
be a spatially small wavepacket which can be treated in the geometric optics approxima- 
tion. The wavepacket 's motion is the same (up to parametrically small corrections) as that 
of a 5-dimensional "particle" traveling in the AdSs-Schwarzschild geometry, and so it fol- 
lows a geodesic whose trajectory is easily calculated in terms of the 4-momentum of the 
excitation. (See section Ull Bl for more detail.) 

It will be important for what follows to remember that the AdS / CFT correspondence is 
really a correspondence between field theory and 10-dimensional string theory. In the strong- 
coupling limit A=oo of the field theory, the general correspondence reduces to one between 
the field theory and the infrared limit of the string theory, which is a supergravity theory. 
The quanta of the supergravity fields correspond to string states that are massless in flat 10- 
dimensional spacetime, such as the graviton. For A=oo, the well-known gravitational dual 
of finite-temperature M = 4 SYM is Type IIB supergravity in an (AdSs-Schwarzschild) x S 5 
background. 



B. What we find 

The classical wave packet falling in fig. H] is a localized, classical excitation of the super- 
gravity fields. For the sake of specificity, consider the case where it is an excitation of the 
background gravitational field. In general, a classical wave is a coherent superposition of the 
corresponding quanta, and so our high-frequency classical wavepacket is a coherent super- 
position of high-frequency gravitons. Since our wavepacket behaves like a particle (in the 
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FIG. 4: Qualitative sketch of the motion though AdSs-Schwarzschild of a wave packet with high 
3- momentum in the x 3 direction. As measured by x°, the particle takes infinitely long to reach 
the horizon. Of special importance is the parametric scale x+ in the fifth dimension, where the 
trajectory turns over and beyond which progress in x 3 rapidly slows to a stop. 



geometric optics approximation appropriate for £ stop <C 4nax), let's follow just one of these 
gravitons. So think of trajectory in fig. H] as that of a single high-momentum graviton with 
a localized wave function. As we will discuss later, interactions between the gravitons that 
make up the wavepacket are very small, and so it is adequate to think about the evolution 
of individual high-momentum gravitons. 

A graviton is really a tiny loop of string whose internal degrees of freedom are in their 
ground state. Because of the gravitational field from the black brane, this closed string will 
feel tidal forces as it falls, which will try to stretch the string in some directions and squeeze 
it in others. As the graviton gets further from the boundary (and so closer to the black 
brane), the tidal forces will increase, and eventually they will become large enough to excite 
the internal string degrees of freedom of the graviton. We will argue that it is the excitation 
of these string degrees of freedom that is responsible for the breakdown of the expansion in 
fig. [3] in the problem region (jl.2p . 

We will find that, in the problematic case ( II ,2p where £ stop <C A _1 / 6 £ max , the tidal forces 
are strong enough to stretch that loop of string to become classically large before the stopping 
distance is reached. This is why stringy corrections cannot be ignored in that case, explaining 
the breakdown of the expansion in fig. [3J (In contrast, the tidal forces are not strong enough 
to excite the graviton's internal degrees of freedom soon enough when £ stop ^> A~ 1/,6 £ max .) 
Though the resulting classical string loop will be large compared to the size of a graviton, we 
should ask how its size compares to the stopping distance 4to P - We will find that the ratio 
of (i) the stretched, classical string's size in the direction of motion x 3 to (ii) the stopping 
distance £ stop is parametrically of order 

•E )string \ — 1/4 ( ^ ^ ^max 

4top \ 4top 

Because of the A _1//4 , this ratio is typically parametrically small for large but finite A, and we 
will argue that the stretching of the graviton into a string (and the accompanying breakdown 
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FIG. 5: The (pink) shaded area represents a narrow region of space-time around the null geodesic of 
fig- HI The AdS5-Schwarzschild metric in this region may be approximated as a pp-wave background 
for the purpose of quantizing a small, falling loop of string that describes a graviton (or other 
particle) in the initial falling wavepacket. 



of the formal expansion in 1/vA in fig. [3]) then has sub-leading impact on A=oo results for 
the stopping distance. But (jl.3]l also includes cases where the stretching of the string may 
play an important role: If one considers a situation where the argument of the logarithm in 
( II. 3p is exponentially large, then the logarithm can be large enough to compensate for the 
factor of A _1//4 . We discuss this situation further in our conclusions, where we make contact 
with folded classical string configurations that were originally used by Gubser et al. jEf to 
model jet stopping. 

Since the tidal forces stretch a quantum string (the graviton) into a larger classical string, 
one may wonder whether or not it is possible to do a real, detailed calculation of the transition 
between the two. Having restricted attention to a single graviton, our problem reduces to 
following the evolution of a single closed string in the AdSs-Schwarzschild background. In 
general, it is not known how to quantize a string in an AdSs-Schwarzschild background. But 
remember that our graviton is localized and so only probes a region of space-time near the 
geodesic depicted in fig. HJ It is enough to consider only a narrow region of the space-time 
that lies near a null geodesic, as depicted in fig. O and so we may treat the full background 
metric in an approximation (known as a Penrose limit) that treats displacements from the 
null geodesic as small. The resulting approximation to the background metric is an example 
of what is known as a pp-wave background, and it is known how to quantize a string in a pp- 
wave background. In particular, it will be possible to calculate the probability distribution 
of the shape of the classical string loop. The methods we use are similar to previous works 
by other authors on the excitation of string modes in scattering processes and/or in pp-wave 
backgrounds 

0-113. 



In the next section, we give a simple, back-of-the-envelope argument that yields the 
result (5x 3 )/£ stop ~ A~ 1//4 but is not precise enough to explain the logarithmic factor in 
( II. 3p . Back-of-the-envelope estimates are sometimes enlightening and sometimes frustrat- 
ingly unconvincing, and the rest of the paper is devoted to a more formal calculation along 
the lines just described. We set up our notational conventions and review A=oo calcula- 
tions in section IIHI Next we review in section [IV] the source of the breakdown of the a' 
expansion for 4to P *C A _1 / 6 £ max , depicted in fig. [3j and use it to motivate why these prob- 
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lems should be overcome by following the quantum evolution of single closed strings in the 
AdSs-Schwarzschild background. In section |V] we take the relevant Penrose limit of the 
AdS5-Schwarzschild metric, crudely depicted in fig. [5j We are then ready to quantize the 
string in section IVT1 and solve for the evolution of the graviton into a classical loop of string. 
We then compute the average size of the resulting classical loops of string, which parametri- 
cally gives fjl.3j) . including the logarithmic factor, for the interesting case £ stop <C A _1 / 6 £ max . 
Some readers may wonder how a quantum treatment of the stretching of gravitons can un- 
derlie the analysis, given that the gravitational/string dual theory is supposed to be classical 
in the limit N c =oo taken in this paper. We discuss this in section IVHt where we also give 
a more detailed justification for treating the gravitons in the wavepacket as independent. 
In section IVIIII (supplemented by Appendix [Bj , we then revisit the Penrose limit used in 
our analysis and verify that it is justified, provided that (I1.3P is small. Finally, we offer our 
conclusions in section IIXI 

II. A BACK-OF-THE-ENVELOPE ESTIMATE 

In this section, we will make a parametric estimate of the amount of tidal stretching of the 
string compared to the size of the stopping distance £ s top- In a later section, we will review 
the A=oo results for how the stopping distance depends on the energy and 4-virtuality of 
our jet source, but here the only thing we will need to know is that the stopping distance 
given by following a null geodesic as in fig. [4] is proportional to a power of the slope dx 3 /dx 5 
of that geodesic where it starts, at the boundary. The more downward-directed one starts 
the trajectory in fig. 0J the less distance it will travel in x 3 before reaching the horizon. 

Now interpret the trajectory of fig. H] as a trajectory for the center of mass of a tiny, 
falling loop of string. Once the string gets far enough from the boundary that the tidal 
forces dominate over the string tension, then the string tension becomes ignorable, and 
different pieces of the string will fall independently along their own geodesies, the string 
stretching accordingly. Imagine plotting two such geodesies, for the two bits of the string 
loop that are most separated. The separation of those geodesies is a measure of the extent 
of the tidally-stretched loop of string as it falls towards the horizon. The proper size of the 
string should start out of order the quantum mechanical size E of the graviton, which is 
roughly set by dimensional analysis in terms of the string tension T as 



where a' = l/2nT is the string slope parameter. 

Very close to the boundary, the tidal forces due to the black hole are negligible, and the 
closed loop of string is in its ground state. We can set up our two geodesies above so that, 
correspondingly, they maintain constant proper separation E graviton near the boundary, where 
the AdSs-Schwarzschild metric approaches a purely AdSs metric. To see how this works, 
imagine making a 4-dimensional boost from (i) the plasma rest frame, in which we create 
an excitation with large 4-momentum = (ou, 0, 0, q^) ~ (E, 0, 0, E) and relatively small 
4-virtuality — q 2 <C E 2 , to (ii) the excitation's initial rest frame, where the 4-momentum is 
instead (y — q 2 , 0, 0, 0). The Lorentz boost factor for this transformation is 




(2.1) 




(2.2) 



S 



stopping distance / slop 




FIG. 6: (a) Parallel geodesies in AdSs in the Lorentz frame where the excitation is at rest in 
3-space. These geodesies maintain a constant proper separation as they fall into the bulk, and this 
separation should be thought of as of order the characteristic size (~ Vo/) of the closed quantum 
string loop describing the graviton (or other massless string mode). The narrow red loops are 
meant to be suggestive of the closed string loop, (b) The same picture boosted to the original 
Lorentz frame, (c) A picture of how those geodesies evolve in AdSs-Schwarzschild rather than 
AdSs. The early-time behavior is the same as (b). [For classical oscillating string solutions, the 
strings depicted in (a) may be thought of as snapshots at moments when the string's proper extent 
in a; 3 is at, say, maximum (or half-maximum or whatever). Such solutions would similarly oscillate 
in the z <C z+ part of (b) but not in the z>z 4 part, where tidal forces dominate over tension.] 



In AdSs, the trajectory in the new frame will drop straight down away from the boundary, 
as depicted by the dashed line in fig. [Bh, 

Now consider the graviton as an extended object with proper size S. The two straight 
solid null lines in fig. [6^ depict the extent of the graviton in AdSs in the excitation's rest 
frame at early times. In Poincare coordinates for pure AdSs, 

R 2 

(ds) 2 = — (dx^dx" + (dzf) , (2.3) 

null geodesies are straight lines. Here R is the radius of the 5-sphere S 5 . We parametrize 
the two solid lines of fig. as 

x 1 = ( 7+ ,0,0,±/VY+,l)2 (2.4) 

with f3 + <C 1 and 7+ = (1 — Z? 2 ) -1 / 2 — 1. Because of the warp factor in the metric, these 
two lines are parallel and maintain constant proper separation 



V Ax 3 £33 Ax 3 = 2/3+7+R ~ 2/3+R (2.5) 

as a function of the rest-frame time. Setting this proper separation to be of order the 
graviton size £ given by ( 12. ip then gives 



_ graviton _ V«_ (2 g) 

+ R R V 1 

The AdS / CFT dictionary relates a' in the string theory with the 't Hooft coupling of Af—4 
SYM by O 

^ = A" 1 / 2 = (£ YM Ag- 1/2 . (2.7) 
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So 

fi+ ~ A- 1 / 4 , (2.8) 

and (12. 4p gives 

x 1 ~ (l,0,0,±A~ 1/4 ,l) ^. (2.9) 

Now boost back to the original plasma frame using ( 12. 2 j) to get the early-time trajectories 
depicted by solid lines in fig. [6b: 

~ ( 7 (1 ± A" 1 / 4 ), 0, 0, 7 (1 ± A- 1 / 4 ), l) z, (2.10) 

where we have used 7 ^> 1 (12.21) . Then 

dx 3 



dz 



^ 7 (1±A- 1/4 ), (2.11) 

initial 

which represents a small relative variation 

A(dx 3 /dz) 



dx 3 /dz 



A~ 1/4 (2.12) 



initial 

in the initial slope dx 3 / dz of the trajectory. As discussed before, the stopping distance (which 
requires a calculation in the full AdSs-Schwarzschild metric) covered by a null geodesic is 
power-law related to this initial slope, and so the difference A£ s top i n how far the two bits 
of string travel also has the same small size ( I2.12p relative to £ s top'- 

Af 

stop - A" 1/4 . (2.13) 



*-stop 

This is just our parametric result (I1.3P quoted in the introduction but without the logarith- 
mic factor. The logarithmic factor requires a more detailed analysis. 



III. SETUP 
A. Notation 

In this paper, we will use Greek letters for 4-dimensional space-time indices (fi, v = 
0, 1,2,3) and upper-case roman letters for indices that run over all 5 dimensions of AdSs- 
Schwarzschild (I, J = 0,1,2,3,5). One form of the metric we use for AdS 5 -Schwarzschild 
is 

(ds) 2 = % [-f{dtf + (dx) 2 + r\dz) 2 ] , (3.1) 

z z 

where z is the coordinate x 5 of the fifth dimension, R is the radius of the 5-sphere (which 
will drop out of final results), and 

z 4 



7=1-3T- 
The boundary is at z=0, and the horizon is at 



-A 



We will not need to worry about the details of regularizing the location of the boundary in 
this work. 
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B. Review of A=oo calculation 



Fig. [T] gave a cartoon picture of how we create our "jets" in the plasma. Readers may find 
a precise description of the field theory problem in any of the previous papers jl, 0, 11, 17 



utilizing this method, but we will not need those details here. Suffice it to say that, in the 
gravity description, the boundary is perturbed in a localized region of space-time in such 
a way as to create an excitation with large 4-momentum ~ (E, 0, 0, E) and a relatively 
small amount of time-like 4-virtuality — q 2 <C E 2 . The response of the system is then tracked 
at late times to see where the excitation comes to a stop. 

The supergravity field that is excited depends on the nature of the source for the jet 
and is the supergravity field dual to the vertex operator in fig. [TJ For example, if we study 
the decay of a slightly off-shell graviton in the 4-dimensional quantum field theory, then 
the relevant supergravity excitation is in the 5-dimensional gravitational field; if we were to 
study the decay of a gauge boson weakly coupled to R charge in the 4-dimensional quantum 
field theory, then the supergravity excitation would be in a corresponding 5-dimensional 
gauge field; and so forth. 

As mentioned in the introduction, there is a nice simplification to the A=oo analysis of 
this problem j§| when £ stop <C £ maiX - On the gravity side, the wavepacket's motion is then the 
same (up to parametrically small corrections) as that of a 5-dimensional "particle" traveling 
in the AdSs-Schwarzschild geometry, and so it follows a geodesic whose trajectory is easily 
calculated in terms of the 4-momentum g M of the excitation. The exact geodesic depends 
on the mass of the 5-dimensional "particle" and so on the mass m of the 5-dimensional 
supergravity field that we have excited, but this mass (if any) may be ignored in the high 
energy limit. 5 

As a result, attention may be restricted to null geodesies, which are given by 

I" n^ v n 
„5\ _ / r— 7 5 9 1" 



Using the metric (13.11) . the A=oo stopping distance is then found to be (9[ 



'Stop — 







* dz M ~ i!® (*-\ V4 — (3 5) 

I 2 + ^i |2 - (4vr)V2 \- q 2 ) 2vrT ' ^ 



An important feature of the integral in ( 13. 5 p is that it is dominated by small values of z, of 
order 

f-q 2 \ l/A f-q 2 \ 1/4 , , 

z± correspond to the parametric scale x\ in fig. HJ where the trajectory turns over and beyond 
which 3-space motion rapidly slows to a stop. The stopping distance is determined by the 
behavior of the trajectory at z ~ z*. 



5 The relevant 5-dimensional supergravity field is the one dual to the operator used to create the "jet" 
excitation in J\f=A SYM. The mass m of the supergravity field is determined by the conformal dimension 
A of that operator, e.g. (Rto) 2 = A(A — d) for scalar operators [18j . where d=A. We emphasize that m 
is a mass in the 5-dimensional supergravity theory and has nothing to do with "mass" of a jet from the 
point of view of the 3+1 dimensional Af=4 SYM quantum held theory. 
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IV. DISCUSSION OF EXPANSION 



We now want to consider corrections to the A=oo results for the stopping distance in the 
case N c = oo. First, we take a moment to review the generic story of l/\/A corrections in 
the AdS/CFT correspondence, which relates 19l 



^ string 4vriV c ™' (41) 
' (i&tNj-W (4.2) 



2ttTR 2 R 2 " wm w ' 

where <? s t r ing is the string loop expansion parameter. The string tension T sets the mass scale 
for massive string excitations, and so a' — > corresponds to taking the scale for massive 
string excitations to infinity. For A=oo, the strongly-coupled 4-dimensional quantum field 
theory is therefore dual to the infrared limit of the 10- dimensional string theory, namely 
supergravity, in the appropriate background. For large but finite A, massive string modes 
are not completely ignorable, and the effective supergravity theory of the massless modes 
gets corrections, in the form of higher-dimensional terms in its action, from integrating 
out the effects of the massive modes. Schematically, the effective supergravity Lagrangian 
becomes 6 

£ ~ R + [a' 3 R 4 + a' b D 4 R 4 + a m D 6 R 4 + • • • ] 

+ [a' 5 D 2 R 5 + a' 6 D 4 R 5 + a /7 D 6 R 5 + ■ ■ • ] + ■ • • , (4.3) 

where we have focused just on the gravitational fields for simplicity. R represents factors 
of the Riemann tensor, and we have not shown numerical coefficients or how the indices 
contract. For N c —oo, there are no loop effects (<7string=0), and accounting for the massive 
string modes in the effective theory is analogous to replacing the effects of the W boson by 
the Fermi 4-point interaction in electroweak theory. So, for example, the R 4 terms in (14. 3 j) 
are calculated from string amplitudes for 2—7-2 graviton scattering and, crudely speaking, 
they correspond to processes which involve intermediate massive string states, as depicted 
schematically in fig. [7J (A more accurate statement would be that they correspond to the full 
string amplitude for graviton-graviton scattering minus the sum of the s, t, and ^-channel 
graviton exchange diagrams that one would calculate in the a'—O supergravity theory.) The 
R 5 terms similarly account for corrections to the 5-point graviton interaction, and so forth. 

In our application, we are interested in the evolution of a high-energy excitation prop- 
agating through the soft AdSs-Schwarzschild background. For simplicity of presentation, 
we will focus on the case where the excitation is in the 5-dimensional gravitational fields, 
though our conclusions will not be sensitive to this assumption. A classical gravitational 
excitation can also be thought of as a coherent configuration of gravitons, and the rele- 
vant string scattering amplitudes are those where two of the external lines are the incoming 
and outgoing high-energy gravitons and the others are the soft background field. So, for 
a 4-point scattering amplitude such as fig. [7J the relevant kinematic limit is that depicted 



6 The precise details of which terms D m R n appear independently in (|4.3|) and which do not will not matter 
to our discussion, but Table I of ref. [20] gives a nice summary of what's currently known at tree level 
(i.e. corresponding to N c =oo). 
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graviton 




graviton 



graviton 



FIG. 7: A picture of massive string mode corrections to graviton-graviton scattering that are 
accounted for by the D m R 4 corrections to the effective supergravity action. 
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FIG. 8: A high-energy graviton, depicted as a string loop, interacting twice with the AdSs- 
Schwarzschild background gravitational field. 



in fig. [SJ With the notation used in that figure, the high-energy limit corresponds to po- 
tentially large s = — (pi + P2)i(pi + P2) 1 but small t = — (pi — Pz)iipi — P3) 1 ■ [Here and 
throughout we may think of the p's as 5-dimensional momenta in AdSs-Schwarzschild rather 
than 10-dimensional momenta in (AdSs-Schwarzschild) x S 5 because in our problem there is 



no interesting dynamics associated with the 5-sphere S 5 .] As discussed in ref. [UJ], the 
D m R 4 terms in (gSD all become equally important in the jet stopping problem when this 
5-dimensional \fs becomes large enough at z ~ z+ to excite massive string modes in fig. |HJ 
The string mass scale is of order 1/a/o 7 , and this condition 

V^d^ £ A= (4-4) 



is shown in ref. [rjj to be the same, in the context of the jet stopping problem, as the 
condition 

which is the problematic case (II. 2p highlighted in our introduction. In this region, massive 
string states in the intermediate state in fig. [H] are kinematically accessible and cannot be 
ignored. 

As the high-energy excitation falls from the boundary to the horizon, as in fig. HI it does 
not just interact with the background field once or twice but does so over and over again, 
as depicted in fig. |9j If the massive string states are kinematically accessible as in ( 14.4)) . 
then they cannot be neglected in any of the internal lines, which means in the effective 
theory language of (14 . 3 1) that all D m R n terms will also become important. This is just what 
happens at the £ stop ~ A~ 1 / 6 £ max point in fig. [31 where all the corrections become the same 
size, corresponding to the threshold ^/S( 5 .di m ) ~ l/va7. Also, note that if massive string 
modes are kinematically accessible for intermediate states, then they are also accessible as 
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soft soft soft soft 

background background background background 



FIG. 9: Like fig. [8] but with many interactions with the soft background field. 



final states. (In our problem, however, there is never really an ultimate "final" asymptotic 
state of the graviton because the graviton falls into the black brane.) 

Interaction with a pure AdSs background does not produce these massive string mode ex- 
citations. On a technical level, the i? 4 interactions really involve the Weyl tensor (the trace- 
less part of the Riemann tensor), which vanishes for AdS5 but not for AdSs-Schwarzschild. 
So it is only the gravitational effect due to the presence of the black brane that contributes to 
massive string mode excitation. As a result, the effects of excited string modes are negligible 
at the boundary and become stronger as one moves away from it (and so closer to the black 
brane). At some distance from the boundary which we will review later, the gravitational 
effects of the black brane become strong enough that (14.41) is satisfied, which is when string 
modes may first be excited. 

From the point of view of an effective theory (14.31) of gravitons, having all the correction 
terms become the same size (or worse), seems like an hopeless disaster for the purpose of 
computations. However, the picture of fig. M suggests a different tack. What is happening 
is that the 10-dimensional gravitons which make up the classical excitation are really tiny 
(quantum) loops of string which are getting their internal string degrees of freedom excited 
as they fall in the background gravitational field. Specifically, internal degrees of freedom of 
a small object are affected by gravitational tidal forces, which try to compress the object in 
some directions and stretch it in others. In any case, consider the fate of a single graviton 
as depicted by fig. M a high-momentum object moves through a soft background field. 
Various authors have previously studied applications of the eikonal approximation to string 
scattering [H, M ■ The upshot, as reviewed below, is that fig. M may be replaced by the 
evolution of a single string quantized in the classical background field. 

One might simply assert that the right thing to do in the eikonal limit is to quantize a 
single string in the classical background field, but there is a very nice paper by D'Appollonio, 
Di Vecchia, Russo, and Veneziano 12| that explicitly checks this in a closely related context. 
As the source of their gravitational field, they take a stack of N coincident Dp-branes 
in otherwise-flat 10-dimensional space-time. They then probe this gravitational field by 
scattering a high momentum, massless closed string (such as a graviton) from it. The 
geometry of the situation is depicted in fig. [TD1 which shows the probe particle moving 
in two of the asymptotically-flat spatial directions perpendicular to the Dp-branes. The 
particle is deflected by the gravitational field of the Dp-branes. (If these were DO-branes in 
4-dimensional space-time, this could be a picture of a particle deflected by the gravitational 
field of the Sun.) Their problem is slightly different from our problem in that their particle 
eventually escapes back to infinity, given the geometry of their setup, but never mind that. 
They further assume that the impact parameter b in fig. [10] is large enough that bE 3> 1. 
This is the assumption that the background field experienced by the particle is soft, its 
momentum components ~ 1/6 small compared to the particle's energy. After discussing the 
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Dp-branes 



FIG. 10: A high energy massless string mode, such as a graviton, deflected by the gravitational 
field sourced by a stack of Dp-branes. The plane of the figure is a plane orthogonal to the Dp-branes. 
(So, for instance, a Dl-brane could be visualized as a line extending out of the page.) 



high energy 
graviton 





high energy 
graviton 



Dp-branes 



FIG. 11: The 2nd-order scattering amplitude for a graviton to elastically scatter from a stack 
of Dp-branes, calculated as a string scattering amplitude with two connections to the Dp-branes. 
This is only a topological picture: As described in fig. [T0| the motion of the graviton is orthogonal 
to the Dp-branes in the problem studied by D'Appollonio et al. 



12 . 



eikonal approximation more generally, they then make the following check. Consider the 
elastic amplitude for the string loop to interact exactly twice with the Dp-branes as it flies 
by and to emerge in the same massless state at the end. They perform this calculation in two 
different ways. One way is to do the full string calculation in the presence of the D-branes, 
as depicted in fig. [TTJ The other way is to simply quantize a single string loop in the classical 
gravitational background caused by the D-branes by taking the Penrose limit of the metric 
near the string trajectory and quantizing the string in the resulting pp-wave background. 
Then they calculate to second order in that background. In the eikonal limit, they verify 
that they get exactly the same result with either method. The calculated probability for the 
string to remain in its massless state drops rapidly below 1 once the kinematic threshold for 
exciting internal string modes is exceeded. 

With this reassurance, we now turn to taking the Penrose limit and quantizing the string 
in our own problem. 

V. THE PENROSE LIMIT 



We begin by taking the Penrose limit to describe a narrow region around the null geod esic 



03.4)) . depicted earlier in fig. |5j 7 For an overview of taking Penrose limits, see Refs. 23M25 



Penrose limits have previously been studied in AdSs-Schwarzschild by Pando Zayas and Sonnenschein 



221 ] . but the null geodesic studied was different. Their geodesic fell straight toward the horizon in AdSs- 
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The null geodesic (13. 4 p can be written as 



dx" 



g" u g v 

OJ 



du, 



where u is an afhne parameter for the geodesic determined by 8 

/g~55 dx 5 



du 



OJ 



(5.1a) 



(5.1b) 



which can be integrated to give u as a function of x 5 . The normalization of u is just 
convention, and the canceling factors of oj = —go in (15. ip have just been chosen to give u 
dimensions of length. In the metric (13. lj) . u is given by 



u 



R 



dz 



z 2 {l- f\q\ 2 /0J 2 y/ 2 



(5.2) 



which is linearly divergent as z — > 0. Our convention will be to define u as running from 
u(z=0) = — oo at the boundary to u(z=oo) = at the black brane singularity, and so 



u{z) 



-R 



dz 



z 2 (l- f\q\ 2 /0J 2 y/ 2 ' 



(5.3) 



The important thing to remember in what follows is that late times correspond to small 
negative values of u (which we will later also call r). Also, though it will be convenient to 
have defined u = to be at the singularity z — oo, all of the physics we discuss will only 
depend on what happens outside of the horizon, z < z^. 

We'll refer to our reference geodesic as x"(x 5 ), which we choose to start at the origin 
x" = on the boundary. Now measure 4-positions in AdSs-Schwarzschild relative to this 
geodesic by defining 

Ax" = x» -x»{x 5 ). (5.4) 
Take q to be in the x 3 direction. Changing coordinates from x 5 and Ax° to u=u(x 5 ) and 



qy. Ax" 



OJ 



-Ax° + — Ax 2 

OJ 



(5.5) 



puts the AdS 5 -Schwarzschild metric (13. ip into the form 



R 

(ds) 2 — 2 du dv H — - 

z 2 



(dx 1 ) 2 + (dx 2 ) 2 + 



(^ 2 -/|<zl 



3\2 



OJ" 



(dAx 



+ 2f^dv dAx 3 

OJ 



f(dv) 



, (5.6) 



Schwarzschild, corresponding to q=0 in our problem rather than \q\ ~ E. Their geodesies also have 
non-trivial motion on the 5-sphere S* 5 . In our application, no dynamical evolution of the S 5 degrees of 
freedom takes place, and the S 5 string degrees of freedom simply remain in a quantum state given by the 
S 5 harmonic of the supergravity field of interest. 
8 This u should not be confused with the coordinate u = (z/z^) 2 used in earlier work by some of the 
authors HQ. 
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where / is now implicitly a function of u. The Penrose limit consists of keeping only the 
terms in the metric that would dominate after a scaling of coordinates 

u -» u, v^j~ 2 v x % -> 7~V (5.7) 

for very large 7. This is analogous to making what would be a very large boost (u — > ju, 
v — > 7~ 1 t> ) in flat space, and so looking at physics close to the light cone, followed by 
rescaling all coordinates by a factor of 7 _1 . For ( 15. 6p . the resulting limit is 



R 2 

(ds) 2 = 2 du dv H — - 



(dx 1 ) 2 + (dx 2 ) 2 + ^ (dAx~ 



(^ ~f\qF' 
,2 



(5.8) 



which is a particular example of a pp-wave metric. When we are done with our analysis, we 
will go back and check that this approximation is justified. The coordinates used in (15. 8p 
are known as Rosen coordinates. 

The metric (I5.8P has the schematic form 

(ds) 2 = 2dudv + ^Kj(tt) (dyi) 2 . (5.9) 

i 

It is useful to normalize the last term by switching to Brinkmann coordinates 

Vi = ViV K i( u ) ( 5 - 10 ) 

and 

v = v-\Y,9u{^V^)y 2 i (5-11) 

i 

to give 

{dsf = 2dudv+ 5^) 2 + (J2 %r ^) {du)2 • (5,12) 

i i » 1 

In our case, using (15. 2\ to rewrite 

^ = ^(W^) 17 X (5-13) 
the metric in Brinkmann coordinates is 

ds 2 pp = 2dudv + (dx 1 ) 2 + (dx 2 ) 2 + (dAx 3 ) 2 + Q(u, x\ x 2 , Ax 3 ) (du) 2 (5.14) 

with 

g(u,x\x 2 ,Ax 3 )=g 1 (u) [(i; 1 ) 2 + (a; 2 ) 2 ] + g 3 (u) (Ax 3 ) 2 , (5.15a) 

d 2 (z~ 1 ) z 3 f'\a\ 2 z 6 \a\ 2 z 6 

ft( „ ) = W „ ) = ^_i = ^L = _ 2 ^_ !!! _ 2 ^, (5,5b) 

= Sj^-m = Af^TM = 4^1. 4^. (5,5c) 

2R 4 w 2 ^R 4 w 2 z 4 R 4 



Here, primes denote derivatives with respect to z, and z=z(u) is implicitly a function of 
m, determined by inverting (15. 3p . The metric ( I5.14p would be flat if not for the Qi, which 
arise from tidal forces. Note that these tidal terms vanish for null geodesies in pure AdS 
(/=1, or equivalently z^— >oo), in agreement with general arguments that the Penrose limit 
of AdS is flat Minkowski space [23] . They also vanish in AdSs-Schwarzschild if q=0 (i.e. if 
the excitation fell straight down in the x 5 direction as a function of time a; ). 
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VI. QUANTIZING OUR FALLING STRING LOOP 



Overview 



We may now follow other authors to quantize a string in a pp-wave background 12l 16 
The bosonic sector is described by a a-model in the pp-wave background metric: 

S = J drday^^idaX^X^guiX), (6.1) 

where 7 is the world-sheet metric and X 1 are the world-sheet fields corresponding to the 
coordinates. For the pp-wave space-time metric (15.141) . this takes the form 

S = So ~aL ; j dTda ^^(daU)(dpU)g(U,AX), (6.2) 

where So is the Minkowski string action. Identifying world-sheet time r with the affine 
parameter u for the pp-wave space-time metric (I5.14p then gives a constraint equation for 
d a V and gives the light-cone gauge Lagrangian 



n u P dn / 

L = — — V ( d T AX i d T AX { - (a'p u )- 2 d a AX i d a AX i + ft(r) AX 1 AX 1 
2 Jo 2n i K 

00 

E [d T AxC d T AP n -ul n {r)AxC AP n ) (6.3) 

i n=— 00 

for the AX\ where p u = v% ^ E and AX 1 = ]T n AX l n (r) e in<T and 



"2 



(We have suppressed the fields corresponding to S 5 coordinates because they will play no 
role in our discussion.) Note that we have chosen a convention where r has units of time 
and a is dimensionless. 9 

Each mode AX l n of the string is a time-dependent harmonic oscillator problem with 
classical frequency u^ n (r). The AX\ and AX 2 modes are tidally compressed as the string 
moves away from the boundary since Q\ = Q2 is negative in (15.15bj) . so that the curvature 
| ^i,n(' r ) of the harmonic oscillator potential increases with time r=u. In contrast, the AX^ 
oscillators are tidally stretched, since Q% is positive in fl5.15c[) . We will focus on these AX 3 
degrees of freedom and so will focus on 

<n{r) = - T — 2 -G,{r). (6.5) 



9 A more conventional normalization might be to define a dimensionless r as r = u/(a'p u ) instead of our 
r = u. One may convert to this convention by everywhere replacing our r by a'p u r. In our convention, 



the world-sheet metric is 

' -(a'p")" 1 
a'p 
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When the string gets far enough from the boundary (i.e. at late enough times r), the Qs 
term in (16.51) becomes dominant and AX^ oscillators become unstable. Physically, this is 
when tidal forces come to dominate over string tension. Using fl5. 15cj) . this instability occurs 



when z > z n , where 



n^ 2 R 2 \ 1/3 . 1/B fnz^ 11 '' 



=A i/e ^ =n V3 (e.e) 

2a' E J \2E J K J 

The instability for the center-of-mass mode n=0 is not particularly interesting: it has nothing 
to do with exciting internal degrees of freedom of the string and just reflects the slight spread 
of the falling wavepacket in figjU due to curvature effects. That is, it reflects physics already 
incorporated in earlier work [9] on our jet stopping problem and has nothing to do with 
gravitons being tidally stretched into classical strings. So, in what follows, we will focus on 
the n^O modes, and the first tidal instability kicks in at z a z\. 

Recall that z*, defined in (I3.6p . characterizes the scale where the x 3 motion of the bulk 
excitation in fig. H] is coming to a stop (no significant motion for z ^> z*). How does the 
instability scale z n above compare to z*? Using (13.51) and (II. ip . the ratio is 



z n n^\V\zHE)V* nV3 4top 



z* z h (-q 2 /E^ A-V64 



(6.7) 



So the tidal instability kicks in before the stopping distance is reached (z\ < z*) precisely 
when we are in the interesting regime £ s top ^ A _1 / 6 £ max ( 11.21) identified in earlier work 
as the case where stringy corrections become important. This is the case we focus on. 
Correspondingly, the modes which become tidally unstable for z < z± are n < n+ with 

Although the instability develops at z = z n , we will see later that the modes n < do 
not have time to stretch significantly until z ~ z*. This fact is suggested by the geodesic 
picture in fig. [6j For z <C 2*, the impact of the black hole on the evolution in fig. [6b is 
negligible, and so the evolution at those times is well approximated by the pure AdS case of 
figs. [6b, and b, for which the proper size of the string remains constant. 

Once a given mode becomes unstable, the quantum mechanics of that mode will be 
somewhat analogous to a standard quantum mechanics thought experiment: What is the 
longest time that an idealized pencil can be balanced on its tip before it falls? Because of 
the uncertainty principle, the pencil cannot be started simultaneously at rest and perfectly 
vertical, and so it must fall. The pencil might be started in a Gaussian wavepacket chosen 
to maximize the average fall time. But it is clear that, once the top of the pencil has 
fallen a macroscopic distance, classical mechanics will suffice to describe its subsequent 
motion: at that time, its position and momentum may, to excellent approximation, be 
considered as simultaneously well defined. For late enough times t\, the pencil's motion for 
t > t\ is approximately classical, with the only effect of its initial wave packet at t=0 being 
to determine a classical probability distribution for the pencil end's position at t\. The 
corresponding momentum at t\ is given (to excellent approximation) by the momentum the 
pencil would have picked up falling classically to that position from vertical. 

In our problem, we will see that, for the string modes n < n* of interest, the analo- 
gous transformation from a quantum description to a probability distribution for classical 
configurations occurs when z ^> z*. 
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B. Solution of the time-dependent harmonic oscillators 



1. Basics 



The distinction between AX 1 (the difference between X % and the reference geodesic) 
and X 1 does not affect the n ^ modes that are our focus. Similarly for the normalized 
coordinates AX 1 . So we will make our notation a little less cumbersome and henceforth 
write AX l n as simply X l n (for 

Each of the real degrees of freedom v^ReX* and \/2\m.X l n in ( 16 .3p have a harmonic 
oscillator Lagrangian of the form 



| mix 2 — oj 2 {t) x 2 ) . 



(6.9) 



with the translation m — )■ p u ~ E and oj 2 {t) — > u 2 n (r). The squared frequency w 2 (r) 
starts at a non-zero value u 2 (— oo) and then changes with time r. The quantum mechanical 
solution to such time- dependent harmonic oscillator problems has a long history. Useful 
explicit formulas for wave functions may be found in ref . [26[ , with applications to strings 

. 10 In our case, the harmonic oscillators all start in 



in pp-wave backgrounds in refs. [13j, [15 



their ground state (the string state describing a graviton) at early times (r — > — oo) and so 
start with Gaussian wave functions. For a time-dependent harmonic oscillator that starts 
as a Gaussian ip(x) oc exp[— x 2 /Aa 2 } at some time r , one may check that the Schrodinger 
equation 



iip 



2m x 



moj 2 (r) x 2 



is solved by 



i/j(x, r) oc 



cxp 



i X(t) 
2X(*) 



mx 



where the complex- valued function x{ T ) satisfies the classical equation of motion 



X 



with initial conditions 



Xi.ro) 



[T)X 



X{To) 



2mcr 2 



(6.10) 



(6.11) 



(6.12) 



(6.13) 
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151,01, and 



10 The case (|6.11j) of interest to us corresponds to the ground state cases (^=0 or n=0) of refs. 

our x here corresponds to their v* , or u* respectively. However, there are various normalization issues 
in the formulas in all of these references, which all differ from us and from each other by overall powers of 
our x* Ix m their expressions for the final wave function, though most differences are probably typographic 
errors, x* Ix i s J us t a time-dependent but x-independent complex phase, and so such differences will not 
matter for the evaluation of an expectation value in one of these states. But, if one wants explicit solutions 
to the time-dependent Schrodinger equation, then the phase needs to be correct. A simple test of any 
such result is to consider the case where w(r) is constant and verify that tl> n (x, t) reproduces the correct 
time dependence exp[— i(n+ |)wr]. Accordingly, the (x/x*Y m (5.37) of ref. 13 1 should be (x/x*)^ + K 



the [v l hf*\ in (2.16) [(2.15) in the preprint] of ref. (l5j should be [v'/v^+i, and the {u/u*) n / 2 in (3.9) 



26] should be (u/u*) 2 + 4 if one wants to keep track of the overall time-dependent phases and have 



of ref. 

results that solve the Schrodinger equation. 
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In our case, where we start in the early-time ground state, that's 

x(to) = 1, x( T o) = iw(-oo) 



(6.14) 



with t — oo. 

It will be useful to have an expression for the corresponding probability distribution 
\ifj(x, r)| 2 for x. From (16. lip , this is just a Gaussian distribution 



Prob(x, r) 



with width 



X r 



2m Im 



(27T)V2 Xrms ( r ) 

X\ [ i\x\ 2 



(6.15) 



(6.16) 



But x*X ~ X*X is a Wronskian of the two solutions % an d X* °f ( EH2J) and so is time- 
independent and may be evaluated at r=r using (I6.14p . giving 



\x(r)\ 



|x(r)| x rms (-oo). 



(6.17) 



Using f )6.5p . x rms (— oo) corresponds to ■sja! '/2n in our application. 

Our remaining task is to solve the classical equation of motion (I6.12p for x- For the case 
of interest X^, using (13.61) . (15.131) . (I5.15cp . (16. 5p . and oj ~ \q\ ~ E 1 , the x equation may be 
put into the form 



fx 

df 2 
dz 

df 



-4(£ 6 - z«) X , 
z 2 (l + z 4 ) l/2 , 



where 



z 

2* 



z 2 R 2 



Z = Zn = = (I 



n\V3 zfc/zi, 



2^ \-Vezt /3 EV 3 ' 



In these variables, the initial conditions (16. 14ft on x are 

x(fo) = h ^(r ) = 2^ 3 . 
dr 

Note from (II. ip and ( 13. 5ft that, parametrically, 



(6.18a) 
(6.18b) 

(6.19) 
(6.20) 
(6.21) 



(6.22) 



"Stop 



(6.23) 
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and so the smallness of is a specific measure of how far we are into the interesting regime 
(11. 2p of 4to P ^ A _1 / 6 £ max . The modes n < n* of interest to us correspond to £„ < 1. 

Using the above equations, one may now check our earlier claim that the string does 
not stretch significantly in proper size at early times z -C {z -C l). 11 But we are more 
interested in what the string does at late times, on which we now focus. 



2. Late-time behavior 



For z>z t (which is z 1), the z equation fl6.18bl) gives 

z = (-3r )- 1/3 , (6.24) 

remembering that our convention is that r is negative and that t(z=oo) = 0. Note from 
(I6.24p that — f is very small at the horizon z = z^, where — f ~ (zh/z+)~ 3 ~ (— q 2 /E 2 ) 3 ^ ~ 
(4to P ^)~ 3 - Throughout this paper we assume that — q 2 <C E 2 and so 4to P T~ l . 

Plugging (16.241) into the x equation (I6.18aj) yields late-time (z ^> 1) solutions x 
(— f) -1 / 3 and x x (~ t) a ^. The dominant solution will be 

X oc (-r)- 1 / 3 . (6.25) 

Though x is a complex-valued function whose purpose is to track the evolution of the 
wavepacket, exactly the same arguments as above give that a classical trajectory would 
have late-time behavior x oc (— f) -1 / 3 . That means that x oc (— f)~ 4 / 3 , and so x and i both 
become large at late times, justifying a classical description at late times. 12 The classical 
relation between the two is determined by a; oc (— f) -1 / 3 to be 

- ~ -— ~ z 3 . (6.26) 

X 37" 

We may extract the proportionality constant in the late-time behavior (I6.25P by solving 
(I6.18P numerically with initial conditions (I6.22p for x, matching the late time behavior of 
the numerical solution to (16.251) . and repeating the calculation for earlier and earlier values 
of r in order to take the Tq — > — oo limit. Our result is that the late-time behavior is 



11 The argument is a distraction and so relegated to this footnote. For z n <C z <C z+ (£ <C z <C 1), one 
may ignore the £ 6 term in (I6.18a[) and treat z as small. In this case the two independent solutions are 
approximately A ~ (£/z)(l + i z 4 ) and B ~ 1 + | z 4 , with d f ~ z 2 (l + ^z 4 )d s . (See Appendix IB 21 for 
a more general discussion of 2 > z„ solutions, for which the ones here are the small z limit.) To match 
to the oscillating solution for x fa z < ( determined by (|6.22p . the superposition of A and B for z ^> £ 
must be x ~ 0(1)A + 0(1)B, where 0(1) denotes coefficients with magnitude of order 1. We then see 
that x is approximately constant over the time period £ z <C 1. 

12 For example, at late times the exponential in the wavepacket (|6.11|) becomes expos'], where S oc x 2 j (— r). 
The WKB condition \d 2 S\ < (d x S) 2 is satisfied as r -> for x cx (-t)" 1 / 3 . We will see shortly that 
the proportionality constant in ()6.25|) is of order 1 for the modes n < n* of interest. If one keeps track 
parametrically of all the proportionality constants in the exponential exp[ii5], one finds more specifically 
that the WKB condition is satisfied when — f <C 1 (i.e. z ^> 1). 
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FIG. 12: The proportionality constant C(£) in (|6,27[) . which determines the late-time width of 
the probability distribution for the amplitude of a string mode. The sloping dashed curve shows 
the large-£ approximation (|6.28a|) . 



with C(£) given by fig. [12j 13 We show in Appendix |A] that the limiting behavior for large £ 
is 

- ^4 f ° r e ^ (6 - 28a) 

shown as a dashed curve in the figure. In the opposite limit of £ <C 1, our numerical results 
approach a constant 

C(£) ~ 0.6428 for £ < 1. (6.28b) 

From ( 16. 5ft . (16. 17ft and (16.271) . and remembering that the analogs of x are 

and \/2ImX^, the late time probability distribution of mode amplitudes X% is given by a 
Gaussian with width 

\X 3 n \ ^-^%-(^-\' 2 for-f«l. (6.29) 
Using ( I6.24p . that may be rewritten as 

/ / \ 1/2 

l*»L - 3 1/3 C(£„) ^ for z » ^,£ n . (6.30) 



13 p or num erical work, it is mildly convenient to eliminate r and express all of the relevant equations solely 
in terms of z, giving z 4 (l + z A )x" + 2z 3 (l + 2z A )\' = -4(£ 6 - z e )x and x'(*o) = 2i£ 3 /^o ( with -> 0) 
and \ X (z)\ -> 3 1 / 3 C(£) z (as z -> oo). 
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As in ( I6.26p . the corresponding momenta in this classical regime are related to the mode 
amplitudes X" 3 by 

d f X 3 n ^^l^^Xl (6.31) 
-3r z* 

Using (15. 8p and (I5.10p . the conversion between the normalized coordinate Ax 3 and the 
displacement Ax 3 from the reference geodesic in Poincare coordinates is 



which is 



Ax 3 ~ — Ax 3 for z > z+. (6.33) 
zR K ' 

So, from (16. 30 p . the amplitudes of the stretched modes in the Poincare coordinate system 



arc 



! A n| rms - d <AW ^ - (2n) l/2 A1/ 4^ 

for fixed r (and so fixed z) in the classical regime. Using (13. 5p and (13. 6p . this may be written 
as 

|X 3 | .3!W^W r i/4, (6 35) 

for 2; > Z*,£ n . 

Note that fixed-r (i.e. fixed-z) slices of the string worldsheet look different than fixed-x° 
slices of the string worldsheet, which is why our depiction of the string at various times x° 
in fig. [6b were not horizontal. 



C. The size of the string at late times 

Parametrically, the average size (I6.35P of each stretched mode X 3 in Poincare coordinates 
is just the A _1//4 4to P that we found in our back-of-the-envelope argument in section [Til But 
we should now look at what happens if we sum up all the modes to get the total average 
size of the string. A convenient measure of the size scale of the string is the average rms 
deviation from the center of the string, 




(6.36) 



where overlines indicate averaging over the string worldsheet position cr and the angle brack- 
ets indicate averaging over the late-time classical probability distribution for each mode 
amplitude. This is given by 




(6.37) 
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(|X^| 2 ) is just the square of what we called |X 3 | rms in (I6.35p . Combining the limiting forms 
( EZH} with (E35]), and recalling from (lOB that £ n = n 1/3 £i, 

3|2 ^/3V3( 87r )V 2 




I^L-(^h^ 4top ) x< r " (|) (6.38) 



where C(0) is given by f l6.28b[) . The sum in (I6.37P is therefore convergent at large n (more 



on that in a moment) and is dominated by a logarithm coming from n = 1 up to n ~ n+ 
At leading order in inverse powers of that logarithm, 

3V3 4 0FC(O) x _ 17 



(5X 3 ) vms ~ I^Uv^i In A^/^topV 7 ^- (6.39) 

1 U/ 



Using (16. 8p . this may be rewritten as 



ii 



(5X 3 ) rms ~ 0.8660 A^ 4 ^ i n W A ~ 1/6 ^* ). (0. 10) 



-stop 



where a parametric expression for the argument of the logarithm is adequate if we are only 
keeping track of the logarithmic term. 

Ignoring the numerical constant in front, (I6.40p is the parametric result (11 .3p that we 
presented in the introduction. One could go on to evaluate the non-logarithmic corrections 
to (I6.40p . but (I6.40p is good enough for our present purposes. We are mostly interested in 
the parametric size of the answer, so that we can determine whether the extent of the string 
in x 3 remains small compared to the stopping distance 4top in the case of large but finite A. 

Qualitatively, what is the origin of the logarithmic factor in (I6.39P and (I6.40p ? Imagine 
doing the same calculation of (SX 3 ) TIDS for the graviton instead of for the stretched, classical 
string. In the ground state, all the modes X % n have an rms size proportional to n -1 / 2 and so 
the sum ^ n {\X l n \ 2 ) is log divergent in the ultraviolet. This does not mean, however, that 
the energy or momentum carried by the graviton is smeared over infinite spatial distances: 15 
in the string's ground state, the bosonic n^O mode contributions to energy and momentum 
are canceled by the fermionic mode contributions, which we have ignored in our analysis. 
Only the bosonic mode amplitudes, however, can become "big" due to tidal stretching, and it 
is the modes whose amplitudes have grown big that we consider when we make the classical 
approximation at late times. Our expression (16.381) is only valid until n gets so large that 



A comment on randomness: |X 3 | rms , determined by (|6.38[) . is the width of a Gaussian distribution for X 3 

I 1/2 



that is centered on X 3 = 0. In contrast, at leading log order, the tr-averaged extent S = (X 3 — X 3 
of the late-time string is always given by the (SX 3 ) rms of (|6.40l) . That's because, when the logarithm 
is large, the calculation of the logarithm in (|6.40p involves a sum over a large number of modes, all 
contributing to the coefficient of the leading log, and so the probability distribution of S 2 is narrowly 
peaked about its probabilistic average (6X 3 ) 2 . ms = (S 2 ) by the central limit theorem. The n=l mode by 
itself is a substantial piece of the first correction beyond leading-log order, and so randomness will enter 
if one moves beyond a leading-log analysis of S 2 . 
15 For some discussion of the unobservability of the UV log divergence and and how the "size" of a low-mass 
string excitation should be interpreted, see refs. 0. 
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the bosonic mode has not been significantly excited (£ n ~ z, and so n ~ n^z/z^) 3 n+ 
when z 3> z±). For yet larger n the modes will be in their ground state and the canceling 
contribution of fermionic modes to physical quantities of interest will come into play. The 
contribution of such high modes is sub-leading and is simply discarded when we approximate 
the string as a classical string at late times. In the resulting treatment of the classical string, 
there is no ambiguity in what the size of the string means, as highlighted by the convergence 
of the sum of (JESS]). 16 

The upshot is that the logarithm in f!6.39j) arises because of (i) the logarithmic UV 
divergence associated with the bosonic modes in the ground state, combined with (ii) the fact 
that the bosonic mode amplitudes with n <C n± all grow by an equal large factor from tidal 
stretching (and no longer cancel against fermionic modes in their physical consequences), 
while those with n ^> n± do not grow significantly in comparison. 

VII. DISCUSSION OF THE SINGLE GRAVITON APPROXIMATION 

We have followed the evolution of a single graviton as it is stretched into a classical string 
loop. We will now take a moment to discuss in more detail the premise that we may follow 
a single graviton. 

First, a single high-momentum string loop might split into two, but such splitting is 1/N C 
suppressed and so may be ignored in the N c =oo limit that we take in this paper. 

The same argument might be given against the possibility of two gravitons merging if 
not for the fact that the localized gravitational wavepacket describing our excitation in the 
bulk is classical, which means it contains a correspondingly large number of gravitons. We 
need a different argument for why interactions between the gravitons that make up the 
classical wavepacket may be ignored. For this, we need to review in slightly more detail the 
formalism of refs. jsl, |9[ for creating the excitation in the first place. In that formalism, the 
initial photon or W boson or whatever of fig. [1] is replaced by a localized, external classical 
field. Specifically, we add a source term to the Lagrangian, 

C + NO(x)e ihx k L (x), (7.1) 

where M is an arbitrarily small source amplitude, 0(x) is a source operator (corresponding 
to the vertex in fig. [TJ, 

~ (E, 0, 0, E) (7.2) 

is the large 4- momentum of the desired excitation, and Al(x) is a slowly varying envelope 
function that localizes the source near the origin in both x 3 and time. The amplitude M 



We point out in passing that this is not an obscure issue specific to relativistic strings. If one quantizes 
small transverse vibrations of an idealized violin string [L = J dz\\px\ + hTx'f\), the same Yln n ^ 1 
divergence arises in the calculation of the mean-square displacement of the string. But, if you took a 
snapshot in time of an actual vibrating violin string (with classically large amplitude), you would not 
have any real-world confusion about what the mean-square displacement of that string meant. You would 
only get confused if you perversely chose to resolve the (idealized) violin string on such small distance 
scales that you could see the quantum uncertainty of very high modes that had not been excited, and 
for such measurements the classical description of the violin string would be inadequate. If you're not 
interested in such high resolution, then the classical result for the average displacement is an excellent 
approximation and is physically meaningful information. 
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is chosen to be small so that we can treat the external source as a small perturbation to 
the strongly-interacting gauge theory, so that the source will never create more than one jet 
with energy E at a time. On an event-by-event basis, the source will usually do nothing at 
all, but on rare occasions (oc |A/"| 2 ), it will create an excitation with energy ~ E. For small 
enough A/", it will essentially never create excitations with energies 2E, 3E, etc. But 
taking M small also means that the bulk excitation created by the source on the boundary 
can be treated in linear response: the self-interactions of the bulk excitation with itself are 
ignorable. For this reason, we may ignore interactions between the gravitons (or other type 
of particles) that make up the bulk excitation. 

Some readers may wonder how we can focus on gravitons, which are quantum mechanical 
objects, when the dual theory for N c =oo is supposed to be a classical theory. As an analogy, 
consider a classical electromagnetic wave with polarization (e x + e y )/\/2. If we choose to, 
we may think of this classical wave as a coherent superposition of photons which are in the 
quantum state (\x) + \y))/y/2, where \x) and \y) are states corresponding to polarization 
in the x and y directions. Now measure the polarization by putting the wave through a 
linear polarizer oriented in the x direction. As we all know, we may view this classically as 
picking out the e x / \/2 component of the wave, or alternatively we may view it in terms of 
photons as saying that each photon has probability 1/2 of being x-polarized. A discussion 
of classical physics in one description is equivalent to a quantum-mechanical discussion of 
probabilities for the behavior of individual quanta in the other description. 

The same reasoning applies to our description of the classical bulk wavepacket in terms 
of individual gravitons. An important conceptual lesson from this is that our wavepacket 
is not a single graviton that evolves into a single classical string loop with probability dis- 
tribution given by a Gaussian with size (16.401) for each mode n. Instead, the wavepacket 
is a large number of gravitons that independently evolve into a large number of classical 
string loops that independently have that probability distribution. Following the analogy 
with photon polarizations further, one could presumably replace the quantum description 
involving probability amplitudes for the degrees of freedom of individual string loops by a 
classical description in terms of a classical string field theory, promoting the Schrodinger 
wave functional t/>[X£] for a single string to a classical wave functional V4^n] that could be 
used to describe all the physics discussed in this paper. We have not pursued this latter 
approach because we thought that the description in terms of gravitons was simpler, more 
intuitive, and more directly related to the previous literature on tidal excitation. 



VIII. CHECKING THE PENROSE LIMIT 



Now that we have our final answer f)6.40p for the size of the classical string that is produced 
by stretching, we should go back and verify that the string is not so big, or so far away from 
the reference geodesic, that the Penrose limit taken in section |V] breaks down. We need to 
check that the dv dAx 3 and (dv) 2 terms in the AdSs-Schwarzschild metric (I5.6p . which were 
dropped in the Penrose limit, are parametrically small compared to the du dv term, 



\dv dAx 3 \ and 

Z 1 0J z 2 



\dvdAx i \ and —r- (dv) 2 < \dudv\, (8.1) 
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for the string motions that we have found. Dividing both sides by \dudv\ and using u ~ \q\, 
we may rewrite these conditions as 



/R 2 



dAx~ 



du 



and 



/R 2 



di > 



du 



(8.2) 



We check these conditions on the string motion in appendix [Bj where we find that the 
condition on \dAx 3 /du\ is the strongest and requires 



A- 1/4 A/rn^< 1 



(8.3) 



in order for our earlier analysis to be valid. Using our result (I6.39p . this condition may be 
written as 



(5X~ 



stop' 



U) 



That is, the Penrose limit only breaks down if one considers the extreme case (to be discussed 
in a moment) where the string becomes as large as the stopping distance itself. 



IX. CONCLUSION 



In our scheme for creating "jets," we have seen different behaviors in the dual theory 
depending on the virtuality (and so the stopping distance) of the jet. For £ st0 p S> A _1 / 6 £ max , 
the gravitons (or other massless string modes) composing the excitation in the gravity de- 
scription remain gravitons until after the excitation has stopped moving in the x 3 direction, 
and there is no difficulty in using the supergravity approximation for the calculation. For 
4top ^ ^ 6 ^max> each graviton is instead stretched into a classical string loop. However, 
provided that 

A -i/4 ln i/ 2 / A-V 6 l m \ ^ (g i) 

\ ^stop / 

the size of that string remains small compared to the stopping distance £ s top- The string 
remains close to its reference geodesic, and so corrections to the A=oo description of the jet 
stopping are parametrically small (if one only attempts to resolve details on size scales large 
compared to the size of the string). However, if instead 

A -l/4 ln l/2 ^- 1/6 W \^ 1; (92) 
\ *-stop / 

then the string loop will stretch out to be parametrically as large as the stopping distance 
itself. Our quantum analysis of the string breaks down in this case (because of the failure 
of the Penrose limit), but we can see what happens qualitatively by tracking what happens 
to the A~ 1/4 log 1/2 < 1 results as we increase the logarithm towards A 1//4 log 1 ' 2 ~ 1. 

In particular, a nice way to visualize what happens is to follow the classical evolution of 
a closed string that initially starts with a proper size £ of order y/a' In n+, which is roughly 
the initial rms size from the modes n < n+ which become classically excited. Increasing 
the logarithm towards A~ 1//4 -\/ln n* ~ 1 is equivalent to increasing E towards ~ R. Fig. [L_3] 
compares examples of such evolution for the cases (a) A -1 / 4 log 1,/2 <C 1 and (b) A" 1 / 4 log 1//2 ~ 
1. (More details of exactly how we initialized our classical string may be found in appendix 
Q ) The interesting feature of fig. [T3b is that, at late times, the string looks like the original 
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FIG. 13: Examples of numerical solutions of the evolution of a falling classical string loop that 
starts near the boundary with proper size (a) £ <C R and (b) £ ~ R. These are snapshots of the 
string at fixed x°. See appendix [C] for details of the initial condition. 




FIG. 14: Schematic pictures of classical folded strings, (a) A closed folded string produced by 
extreme tidal stretching of a graviton in our method of generating "jets" in the case 

ofA-^logVs > 

1 (|9.2p . and (b) the infinite, folded open string considered by Gubser et al. [jj. In the latter case, 
the trailing string continues to get closer and closer to the horizon as x 3 — > — oo. 



picture advocated by Gubser et al. |5| of gluon jets as dual to the evolution of a trailing, 
folded classical string falling in AdS 5 -Schwarzschild. Our string is a folded closed string, as 
depicted in the cartoon of fig. IT4k . whereas the one studied by Gubser et al. was a folded 
infinite open string, as depicted by fig. H4~b. 17 However, the left end of the string in these 
figures, which is very close to the horizon, does not play a significant role in the effect on 
the boundary theory, and so the physics of these two situations is approximately the same. 
Historically, the original motivation of our own method for posing "jet" stopping problems 



17 More precisely, Gubser et al. first considered a folded open string that stretched out from beyond the 
horizon, as in fig. 1 of ref. [5]. But in actual calculations, they focused on the trailing infinite folded string, 
as in fig. 2 of that reference. 



29 



outlined in the introduction to this paper and motivated by fig. [TJ was to give a precise 
field theory problem in A/=4 SYM that could be solved, beginning to end, using gauge- 
gravity duality. It has not previously been know how to precisely set up a problem in J\f=4 
SYM that corresponds to earlier studies of jets [5-7] that made use of classical strings in 
the gravity description. It is interesting to now make contact between our approach and 
Gubser et al.'s classical string approach, in the particular limit (19. 2p . which can be roughly 
rewritten as 

T" 1 « 4to P < T-^ f^j V3 e-°^ l2 \ (9.3) 

where the T _1 <C 4top is thrown in to emphasize that we've always been assuming — q 2 <C E 2 
and so T _1 <C £ stop throughout, and where O means "of order." Alternatively, in terms of 
the virtuality — q 2 of the source of our "jet," (19.31) is 

E 2 » -q 2 > r 4/3£2/3 A 2/3 e +0(AV 2 )_ (g 4) 

This window of stopping lengths only appears once the jet energy is large enough that 

£»7VAe + °( Al/2 ). (9.5) 

Even though there is a region of overlap (19. 31) of our results with strings that look similar 
to those of Gubser et al., there are still important differences once we get out of this range. 
Gubser et al. found a maximum stopping distance of order T~ 4 / 3 (£^/\/A) 1//3 , as do other 
methods that also model excitations with semi-infinite classical strings in the gravity dual 
0]. In contrast, the types of excitations that we create, through processes like fig. [TJ have a 
parametrically larger maximum stopping distance of order £ max ~ y- 4 / 3 ^ 1 / 3 . 

In this paper, we have studied the case of large but finite A, but we have kept N c =oo. 
It is important to mention, following Gubser et al. [5], that the case of finite N c may be 
qualitatively different. For finite N c , the folded classical string may break, creating many 
string loops. It would be interesting to understand whether or not such breaking would 
impact the description of jet stopping in Af—4 SYM. 

Our basic argument in this paper has been that the breakdown of the 1/a/A expansion 



seen in the earlier supergravity calculation of ref. 11] is explained by the tidal stretching of 
gravitons into classical string loops. Logically, however, to show that tidal stretching was the 
only source of difficulty, we should now do a systematic analysis of sub-leading corrections to 
the string evolution that we have presented here and show that those corrections are indeed 
parametrically small. We expect that they will be, but we leave that for someone to study 
another day. 

Finally, we mention that an alternative to our method for creating "jets" is to produce 
gluon jets as synchrotron radiation from heavy quarks that are forced into circular motion. 
The corresponding jet stopping problem has been investigated for strongly-coupled A/"=4 
SYM by Chesler, Ho, and Rajagopal [lOj. We suspect that problems similar to those treated 
here could be set up in that context as well, with appropriate choices of the heavy quark 
velocity and synchrotron radius for a given temperature T. It would be interesting to 
investigate. 
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Appendix A: Large £ behavior of C(£) 



For large £, the z 6 term in the differential equation (16.18aj) for £ can be ignored until z ^> 1. 
At that point, however, we may use the simple large- z result (I6.24p for z. Substituting this 
into (I6.18aj) gives 



fx 

df 2 



-4 



e- 



-3f 



X, 



whose solution is 



X 



The late-time behavior r — > is 



Ixl 



r(I) 



7rV2^(-f)l/3 



(Al) 
(A2) 

(A3) 



Appendix B: Checking the Penrose limit: details 

In this appendix, we will check whether the evolution of our strings satisfy the conditions 
( 18. 2 p required for taking the Penrose limit. It is illuminating to do this in two different ways. 
First, we will make a rough argument based on following divergent geodesies, similar to the 
style of argument that we used in section [III Then we will give alternative (but physically 
equivalent) arguments based on string-based results and formulas. 



1. Tracking diverging geodesies 

As in section (JTTJ) , let us characterize the string by following null geodesies that roughly 
trace different bits of string and which deviate slightly from our reference geodesic. As we'll 
discuss again later, this approximation amounts to ignoring the tension in the string. From 
the null geodesic formula (13 ,4p and the metric (I3.ip . the x 3 coordinate for such geodesies is 
given by 

dx 3 q 3 

where 



Bl 



^ = ^ = (-1,4)- (B2) 
u 

Remembering that Ax^ = x^ — x^(z) is the deviation relative to the reference geodesic, we 
have 

dAx 3 q 3 q 3 



Expand to first order in Aq 3 = q 3 — q 3 : 

dAx 3 Aq 3 



dz (1 - /g 2 )3/2 
Then using ( 15 .131) (and defining u with respect to the reference geodesic x), 

fR 2 dAx 3 f Aq 3 
z 2 du 1 — fq 2 



(B3) 



(B4) 



(B5) 



31 



Since 1 — fq 2 ~ (z* + z 4 )/z^ the combination f ]B5j) is largest for z < z*, and the dAx 3 /du 
condition in (18. 2 j) requires 

Ag 3 « 4 (B6) 



for the Penrose limit. Use (13. 5 p to relate this to the stopping distance: 



f-stop ~ — ~ z h 



E 2 



2 \ 1/4 



2h 



^-gs) 1 / 4 ' 



(B7) 



so that 



At 



z h Aq 3 Aq 3 



stop 



stop 



Aq : 



3 ^-stop ^ ■ 



;i - g 3 ) 5/4 i - h 

Combining (1B6[) and (1B8[) gives the condition 

A£ s top ^stop 

quoted in ( Eg| . 

Now turn to the condition on dv/du in (18. 2p . The definition (15. 5 p of t> gives 

dv = d(Ax M ) = -d(Ax°) + g 3 d(Ai 3 ), 
and so we need a formula for d(Ax°). The analog of (1B1I) is 

cfe -f^qo 



with expansion 



dAx c 



ft Ag 3 



(B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



dz (1 - /<f) 3 / 2 

Combining (IB4|) . (IBlOp . and (IB12|) . gives dv/du ~ 0. We therefore have to go back and 
make our expansions to second-order in Aq. The result is 



du ^ (Aq ± ) 2 _ (Aft) 2 

dz 



2(1 - M!) 1/2 2(1 - fqlfl 2 ' 



(B13) 



and so 



/R 2 



dv 



du 



2(1 - fql 



The corresponding condition on dv/du in (I8.2p is then 

2 /(^3) 2 



/(A<? ± ) 2 and 



< 1. 



(B14) 



(B15) 



The first condition is strongest for z <^ z^ and the second for z < £*, giving 



|Agj.| and |Ag 3 | -| < 1. 

zi 



(B16) 
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Using (1B8I) . the condition involving Ag~3 becomes 



AP <z- P — 

^^^stop *-stor> o • 



stop 2 



(B17) 



Since <C 2^, this is weaker than the previous condition flB9j) . 

Now turn to the other condition, |Aqr_i_| <C 1 in (IB16j) . To estimate |Aq_i_|, return to the 
arguments of section [Til but now, in the rest frame, include an initial proper displacement 
of the two geodesies in a+ of the same parametric size as the initial proper displacement in 
a; 3 . Following through the argument, one finds 



x 1 ~ 



7(l + /9/9+) ) /3 ± , 7 (/9 + /9+) ) l )z 



with (3± ~ (3 + . Then 



Aq 3 = A— = ± dX " /dZ ~ a 7(g + gjj 



ft 



go 



rfx /^ 7(1 + 0/9+) 7 2 



and 



so that 



A<?j 



. Oi . dx ± /dz 
A— = A - — ~ A 

go dx°/dz 

Aq | / E 2 



/3j 



^/3± 

7(1 + ^+) 7 



|A<? 3 



7 



-g 2 zl 



So, using ( IB8 



|A<?j 



I A?; 



3 "J 



stop ^* 



-stop 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



The condition |Aqj_| <C 1 is therefore the same as the previous condition ( 1B17I) and so is 
also weaker than (IB9I). 



2. String-based formulas 



Another way to get to the same result is to start from formulas based on our string 
analysis in section I VI I We will also need the constraint equation for V generated by ( 16. 2 p in 
t = U gauge, which is 



IE 



(d T AX'Y + ——(d„l\X-)' + ft (AX 
(a 'p u y 



i\2 



Using the definition of (15.111) of -0, the corresponding constraint for V is 



d V = -1 V 

2 ^ "T 1 nr) ( a 'p u ) 2 

where Ki = /c»(r) is the notation introduced in eq. (15.91) with 



a it * \ 2 i 



«1 = ^2 



«3 



R 2 (c 2 -/|q| 2 ) R 2 (zj + z*) 



(B23) 



(B24) 



(B25) 
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The relationship ( 15.1 Op between AX 1 and AX 1 is 

AX 1 = ^IAX\ (B26) 

Alternatively, ( IB24I) is just the constraint equation T TT = that one would derive in Rosen 
coordinates ( 15. 8p . 

Here and throughout, we will take results derived in the Penrose limit but then use them 
to check, after the fact, whether the Penrose limit was valid. We will focus on the evolution 
of the string modes after they become unstable (z ^> z n ). The string equation of motion 
from the Lagrangian (I6.3P is 

3 t 2 aj£; = -u4(t)ax;, (B27) 

where u 2 n {r) = n 2 (a'p u )~ 2 — Qiij). At late times z 3> z n , the tidal terms Qi dominate over 
the string tension terms n 2 (a'p u )~ 2 , and the equation of motion becomes 

dlAXi-QiAXl (B28) 

which it is convenient to rewrite as 

d 2 AK - Q ^§- AX;. (B29) 



Two independent solutions for AX in the tensionless approximation ( 1B29|) are 



A\r) oc yfc{T), B\r) oc v^F) / —rr; • (B30) 



dr 



We will take A % n and B l n to be normalized to be ~ 1 at z ~ z n . The late-time (z 3> z n ) 
string solutions will then be given by some probabilistic superposition 

AP n {r) = <Al{T) + U n B\{r) (B31) 

where a\ and b l n are of order the size scale a/ a'/n of the proper size of the ground state of 
the n-th mode, which characterizes the string for z < z n . 
From (I5.13P or ( I6.18bp . the relationship between r and z is 

-^,W°°W^~^#' 2<<2 * ; (B32) 



z \ z * + z ) | , z » z. 

Given the expressions ( 1B25I) for Kj, the resulting behavior of the solution ( 1B30h is then 

A\ = Al ~ % (B33) 



z 

= ^ ~ T (B34) 

2 > z*. 

^nH z (B35) 
5n ~ <!, (B36) 

2 > Z*. 
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Note that B % n dominates over A l n for z ^> z n . For z 3> z±, we have B 3 oc A 3 oc z, and 
this behavior corresponds to the dominant (— f) _1//3 late-time behavior discussed in section 
IVI B 21 At late times, B 3 differs from a multiple of A 3 by corrections of absolute size of 
order z^/z 4 , which is the sub-leading (— f) 4//3 late-time solution discussed in the main text 
just before ( I6.25p . A simple way to see this is to rewrite the definition of B % in (IB30P as 

^-^>L^-y (B37 » 

where the first term is proportional to A i (r) and the second is (at late times) the sub- leading 
solution. Similarly, at late times (z 3> z±), B 1 differs from a multiple of A 1 by sub- leading 
corrections of order z 2 /z 2 . 



a. Condition on d(Ax 3 )/du 



Now let us investigate the Penrose limit condition (18.21) on d(Ax 3 )/du, which can be 
rewritten as 



/R 2 



Or 



AX 3 

\/«3 



<C 1. 



(B38) 



The A oc ^J~k solutions do not give any contribution to the r derivative above, and so 

Bl 



b 3 d 



/K 3 



a' zt/Rz* _ A" 1 / 4 z 2 zi 



n k 3 



yfo R 2 ( Z 4 + Z A ) ' 



(B39) 



The condition (1B38|) has to be satisfied for every point a on the string. But we can get a 
sufficient parametric condition by taking the rms average over the modes n < n± which are 
excited: 

/R 2 



which then gives 



A-^S/ln 



ru 



fzt 



zt + z 4 



< 1, 



< 1. 



(B40) 



(B41) 



This condition is strongest for z < z± and then gives the condition (18.31) quoted in the main 
text. 

To relate this to the earlier geodesic analysis, recall that AX 1 / = AX 1 and so note 
that (IB39D gave 



\d T AX 3 \ ~\-v*J^T / z * - 



Comparing to (1B5|) then identifies Ag3 as 



r 4 4 A n 

a it ^* \-l/4 f\ ^* *- M 'stop 

h *-stop 



(B42) 



(B43) 



consistent with (IB8I). 
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b. Condition on dv/du 



Now turn to the condition on dv/du. The first term in the constraint formula (IB24I) . 
which is what remains if one ignores the tension terms in that equation, simply reproduce 
the conditions (IB15P found in the earlier analysis based on diverging geodesies. To see this, 
consider the i = 3 contribution of the first term in (IB24[) to the dv/du constraint in 



/R 2 



k 3 (d T 



AI 3 \2 



'«3 



< 1. 



(B44) 



Using flB25l) and flB39l) . the contribut ion from mode n to the left-hand side of (1B44I) is of 
order 

JZ * (B45) 



Summing over the excited modes gives 

\- l/2 ln(n 



f4 



<C 1. 



(B46) 



This is the same as the second condition in flB15j) . with the identification (IB43[) noted before. 
Similarly, using (!B25l) and pfflj) . 



/R 2 



AX- 



< 1 



generates 



A- 1 / 2 lnK)^«l, 



(B47) 



(B48) 



which is the same as the first condition in (1B15I) . using the identifications of (IB22I) and 

dHSD. 

Finally, there are the tension terms in the constraint formula flB24j) for dv/du. One may 
check that these also do not generate any constraints stronger than (IB9j) . 



Appendix C: Simulation details for fig. 1131 

Before we explain the initial conditions used to simulate the string evolution in fig. [131 
we will take a short detour and explain how to generate a string evolution as depicted in fig. 
Ek, 

First, let us notice that near the boundary, the metric is approximately AdS. So the 
string evolves at very early times in a geometry which is AdS. Moreover, in the rest frame 
of the string considered in fig. [5k., if the string has a sufficiently large energy, the string 
excitations probe a geometry which is the Penrose limit of an in-falling massless excitation 
in AdS. In other words, these excitations probe flat space. To show this explicitly, we start 
by expanding around the null geodesic x M = (z,0). Defining, as in equations (15. 31 - 15. 5p . 
Ax^ = — x^(z),v = —Ax°,u = —R 2 /z, yields 

ds=2dudv + > ^ —- . (CI) 

z 2 z l 
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Taking the Penrose limit by scaling the coordinates as in (15. 7p brings the metric into the 
Rosen form 

ds 2 = 2dudv+ Y (rfA f )2 , (C2) 

j=l,2,3 

which is just (15. 8p with — > oo (i.e. / — > 1). While this is not immediately recognizable 
as a flat space metric, a change of coordinates to bring the metric into Brinkmann form, 
x 1 = -uAx i /R,v = v - ^ i (Ax i ) 2 /(2w), yields 

ds 2 = 2dudv + ( rf ^) 2 ' ( C3 ) 

1=1,2,3 

which is the z^ — > oo limit of (I5.14p . In this geometry we consider the evolution of an 
initial string configuration with the string shrunk to a point and with d T X°(0, a) = A, 

d T X 3 (0,a) = B cos a, d T Z(0, a) = yj ' (d T X Q ) 2 - (d T X 3 ) 2 , where we have defined X° and Z 

by it = (X° + Z)/2 and v = X° — Z. The other string coordinates remain zero throughout 
this section. Since the string is point-like at the initial time, the stress-tensor constraint 
T rcr = is trivially satisfied, and the second stress tensor constraint T TT = was used to 
solve for the initial condition on d T Z. If we choose A and B such that B C the string 
has an approximate analytic solution of the form X° = At, X 3 = B cos a sin r, Z w At, or, 
in the Brinkmann coordinates, U ~ At, V ~ 0, X 3 = B cos o sin r. 

Next we switch to Rosen coordinates, and we require that the string is again point-like 
at the initial time r = r c : AX°(T c ,a) = 0, AX 3 (r c ,cr) = and Z{r c ,a) = z c , where we 
introduced a boundary regulator z c w 0. In order to satisfy this set of initial conditions 
we only have to perform a shift in U in the previous Brinkmann solution, and use the 
mapping from Brinkmann to Rosen coordinates: U = —R 2 /Z w A(t — r c ) — R 2 /z c , V ~ 0, 
X 3 = B(Z/R) coscrsin(r — r c ). This is a solution which oscillates, as depicted in fig. [TBI for 
small z. The envelope of these oscillations grows linearly with z and roughly corresponds to 

fig.Ek 

The next step is to perform a very large boost along x 3 starting from the rest frame we 
have been discussing so far. This will give us the appropriate initial conditions used to gener- 
ate fig. [T3R Again, the string starts as a point-like object very close to the boundary: Z{r = 
0,a) = z c = 0.01, X°(T,a) = 0, X 3 (r,a) = 0. We take z h = 1. The initial (worldsheet) 
time r = derivatives have the following expressions: <9 r X°(0, a) = 1300 + 0.02 x 1299 cos a, 
d T X 3 (0,(r) = 1299 + 0.02 x 1300 cos a, d T Z(0,a) = ^(d T X°) 2 f 2 (z c ) - (d T X 3 ) 2 f(z c ), where 
f(z c ) = 1 — (z c / Zh) ■ Here 1300 is the energy of the boosted string and 1299 is the momentum 
along the x 3 direction in the appropriate units. 

Earlier in this appendix, we discussed Penrose limits of AdS in order to motivate our 
choice of initial conditions. However, we now simulate the evolution of our classical string 
from these initial conditions in the full AdSs-Schwarzschild geometry, with no such approx- 
imation, to obtain fig. IT3"a . In order to generate a string evolution to sufficiently late times, 
we followed 0, S] and introduced an appropriate stretching function T, 




with m = —0.07, n — 1, p — 4. Then redefine the worldsheet time r so that the worldsheet 
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FIG. 15: An (approximate) AdS string solution for small z showing the type of string we create 
with our initial conditions, viewed in the 4-dimensional string rest frame. The snapshots of the 
string are equally spaced in r. z c could be any point of the oscillation where the length of the 
string vanishes. The straight line geodesies of fig. [6k, may be roughly thought of as corresponding 
to the dashed lines above, which depict the envelope of the oscillation. 



metric is 



-r o 
o r- 1 



(C5) 



The numerical evolution of the classical string was carried out using Mathematica. 

Lastly, to generate fig. [T3b . the initial conditions were modified as follows. The string 
begins its evolution once more as a point-like object close to the boundary Z(t = 0, a) = 
z c = 0.01, X°(t, a) = 0,X 3 (t, a) = 0. But the initial time r = derivatives are now 
<9 T X°(0,a) = 1300 + 0.6 x 1299 cos a, d T X 3 (0, a) = 1299 + 0.6 x 1300 cos a, and, as before, 
the remaining derivative is obtained from the constraint T TT = 0. 

The main difference between the two initial conditions is that we have increased the 
small parameter, which is a placeholder for a/ (ct'/R 2 ) Inn* from 0.02 in the previous nu- 
merical simulation to 0.6 here. [As a side comment, this parameter cannot be increased to 
arbitrarily large values or else the argument of the square root in d T Z becomes negative.] 
The stretching function used here is of the same form as in ( 1C4I) . but with different expo- 
nents m = —0.08, n = l,p = 4.5. This string quickly stretches in z. Once one of the folding 
points reaches the horizon, it remains frozen at that particular value of x 3 . The net result 
is that we have now generated a string which is folded back, with one folding point at the 
horizon and trailing behind the rest of the string, while the other folding point is still close 
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to the boundary. This string is moving with some momentum in the x 3 direction. 
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